arXiv:l506.02237v 1 [nucl-th] 7Jun2015 
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Background: Realistic nucleon-nucleon interactions induce short-range correlations in nuclei. To solve the many-body prob¬ 
lem unitary transformations like the similarity renormalization group (SRG) are often used to soften the interactions. 

Purpose: Two-body densities can be used to illustrate how the SRG eliminates short-range correlations in the wave function. 
The short-range information can however be recovered by transforming the density operators. 

Method: The many-body problem is solved for 4 He in the no core shell model (NCSM) with SRG transformed AV8’ and 
chiral N3LO interactions. The NCSM wave functions are used to calculate two-body densities with bare and SRG transformed 
density operators in two-body approximation. 

Results: The two-body momentum distributions for AV8’ and N3LO have similar high-momentum components up to relative 
momenta of about 2.5 fm 1 , dominated by tensor correlations, but differ in their behavior at higher relative momenta. The 
contributions of many-body correlations are small for pairs with vanishing pair momentum but not negligible for the momentum 
distributions integrated over all pair momenta. Many-body correlations are induced by the strong tensor force and lead to a 
reshuffling of pairs between different spin-isospin channels. 

Conclusions: When using the SRG it is essential to use transformed operators for observables sensitive to short-range physics. 
Back-to-back pairs with vanishing pair momentum are the best tool to study short-range correlations. 

PACS numbers: 21.60.De, 21.30.Fe, 05.10.Cc, 25.30.-c 


I. INTRODUCTION 

Realistic nucleon-nucleon (NN) interactions are fitted to 
NN scattering data up to the pion production threshold. 
Therefore their properties at short distances r and also their 
off-shell behavior is not completely constrained. All realistic 
interactions include pions to describe the long- and medium- 
range parts of the potential. The short-range part is param¬ 
eterized phenomenologically [1], by the exchange of heavy 
mesons [2], or in chiral effective field theory by (regularized) 
contact-terms [3, 4]. At short distances r < 1.5 fm the differ¬ 
ent interactions induce typical short-range correlations due to 
short-range repulsion and the tensor force which is reflected 
in a depopulation of one-body momentum distributions below 
the Fermi momentum and an enhancement at high momenta 
when compared to single-particle mean-field occupation prob¬ 
abilities [5-10]. 

However short-range correlations are two-body correlations 
and therefore two-body densities in coordinate and momen¬ 
tum space provide the best tool to study these correlations. 
Experimentally short-range correlations have been studied in 
inclusive, semi-inclusive and triple-coincidence reactions, see 
the reviews [11, 12] and references therein. The most direct 
information about two-body correlations can be obtained in 
triple-coincidence experiments where one knocks out a nu¬ 
cleon pair with protons [13] or electrons [14-17] at high- 
momentum transfer. In these experiments one found a domi¬ 
nance of pn over pp pairs at high relative momenta, that indi¬ 
cated the importance of short-range tensor correlations. Ten¬ 
sor correlations also appear to play a major role in ip, d ) reac¬ 
tions at high proton energies [18, 19]. 
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Recent theoretical studies of two-body momentum distri¬ 
butions can explain these observations. In a simple picture 
only close nucleons found in relative S -wave pairs are affected 
by short-range correlations and in a first approximation high- 
momentum components are generated by pairs in deuteron- 
like configurations [20], The two-body momentum distribu¬ 
tions at high momenta can also be connected to the nuclear 
contacts [21]. Detailed few- and many-body calculations al¬ 
most exclusively use the Argonne i’| X (AV18) or Argonne v' g 
(AV8’) interactions [22-26]. Compared to the Argonne inter¬ 
actions, interactions derived in chiral effective field theory [3] 
are regularized with a relatively low momentum cut-off and 
one might expect noticeable differences for the short-range 
correlations. 

From the perspective of nuclear many-body calculations 
short-range correlations are not a desired feature but pose a 
severe problem. One way to address this problem is to include 
the short-range correlations explicitly, as in the correlated ba¬ 
sis function theory [27] or in variational Monte Carlo (VMC) 
calculations [22-24]. Another approach is to use soft phase- 
shift equivalent effective interactions that are obtained from 
the bare interactions by means of unitary transformations like 
Vi ow -k [28], the unitary correlation operator method (UCOM) 
[8, 29, 30] or the similarity renormalization group (SRG) [30- 
32]. In real life the unitary transformations are performed in 
an «-body approximation. This usually means in two-body, 
or for state of the art SRG calculations, in three-body approx¬ 
imation [33-35]. Even if we only have a two-body interac¬ 
tion at the beginning the unitary transformation will induce 
three- and higher-body terms. If the unitary transformation 
acts mainly at short distances and the density of the nuclear 
system is low enough, so that the probability to find three nu¬ 
cleons simultaneously close together is small, transformations 
on the two-body level will be a good approximation. The re- 
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maining dependence of observables on the unitary transforma¬ 
tion can be used to analyze the nature of missing many-body 
terms. 

In a consistent calculation it is not enough to transform the 
Hamiltonian, all operators have to be transformed. Within the 
SRG approach the transformation of long-range operators like 
radius or electromagnetic transition operators has been stud¬ 
ied in [36] in two- and three-body approximation. As one 
might expect the effect on these long-range operators is not 
very large. On the other hand we expect large effects for short- 
range or high-momentum observables. Within the UCOM ap¬ 
proach we investigated these effects in two-body approxima¬ 
tion and with simple trial wave functions for the one-body 
momentum distribution [8] and for two-body coordinate and 
momentum space distributions [37]. The SRG operator evo¬ 
lution for the deuteron was studied extensively in [38] and 
analyzed in terms of a factorization for high-momentum ob¬ 
servables. These ideas were extended for Fermi gases in [39]. 

In this paper we investigate two-body densities for the 
ground state of 4 He in coordinate and momentum space for 
SRG transformed AV8’ [1] and chiral N3LO [3] interactions 
using the no core shell model (NCSM). For the bare AV8’ in¬ 
teraction the NCSM results are not converged and we use in 
this case the method of correlated Gaussians [37, 40, 41], We 
will discuss the two-body densities obtained with both, bare 
and transformed density operators. The SRG transformation 
depends on the flow parameter that controls the softness of the 
transformed interaction. We perform the SRG transformation 
in two-body approximation and use the flow-dependence of 
the results to test the quality of the two-body approximation. 
We will show that for pairs with vanishing pair momentum the 
contributions of three-body correlations are negligible. 

After describing in Sec. II the used methods, two-body den¬ 
sities in coordinate and momentum space are presented in 
Sec. Ill where we also discuss the role of many-body corre¬ 
lations for different observables. Summary and conclusions 
follow in Sec. IV. Technical details about the calculation of 
translational invariant two-body densities in the NCSM frame¬ 
work are presented in the Appendix. 


II. METHOD 


The SRG flow equation for the Hamiltonian H a = T int + V a 
and the corresponding transformation matrix U a are given by 


dH a j.— i dU a — 

—— = [q a ,H a \, —— = —U a r\ a , 
da da 


( 1 ) 


with the generator r/ a taken to be the commutator of the intrin¬ 
sic kinetic energy and the evolved Hamiltonian: 


T/a = (2 fl) 2 [Tin,, H a ], Tint = T - T cm . (2) 


Please note that \J a used here and in Ref. [30] corresponds to 
Jj\ in Refs. [31, 32], 

In this paper the Hamiltonian is transformed in two-body 
approximation: 

H a = f%HU a = F int + VM + ... + Vf 1 * fin, + V f ] . (3) 


The evolution (1) can therefore be performed for the relative 
motion in two-body space. The flow equations are solved on 
a momentum space grid with relative momenta going up to 
k max = 15 fm 1 for the AV8’ interaction. In the end the mo¬ 
mentum space matrix elements are evaluated in the harmonic 
oscillator basis to be used in the NCSM. We will present re¬ 
sults for flow parameters of a — 0.01 fm 4 , a = 0.04 fm 4 (a typ¬ 
ical value used in many-body calculations), and a = 0.20 fm 4 
corresponding to a very soft effective Hamiltonian //„ . 

The many-body problem for 4 He is then solved with the 
SRG transformed two-body Hamiltonian H a using the shell 
model code Antoine [42]: 

H n \V (t ) = £ a |'F a > . (4) 

The results for the ground state energy of 4 He and the two- 
body distributions in coordinate and momentum space dis¬ 
cussed in this paper are well converged within the model space 
CVmax = 16, oscillator parameter hQ = 36MeV). The bare 
AV8’ interaction (a = 0) can however not be converged in 
this space and we use here the results obtained with the corre¬ 
lated Gaussian method [37, 40, 41], 

All two-body information contained in the many-body 
states | v P ff ) can be expressed in terms of the two-body den¬ 
sity matrix that is obtained by integrating over all coordinates 
besides the pair position and the relative coordinates of the 
pair. The two-body density as calculated in the NCSM is not 
translationally invariant due to the localization of the NCSM 
wave function in the origin of the coordinate system. How¬ 
ever the translational invariant two-body density can be ob¬ 
tained from the two-body density in the laboratory system by 
a linear transformation as discussed in Appendix D. In the har¬ 
monic oscillator basis we express the two-body density matrix 
P a qQ q'Q' two_ t ,oc ly basis states 

\qQ) = \q)®\Q) 

= |/? re i; nlm, SMs , TMj) ® |/? pa j r ; NLM) . 

Here q summarizes the harmonic oscillator quantum numbers. 
nlm stand for the relative coordinates, S Ms and TMj give the 
total spin and isospin of the pair, respectively, and Q sum¬ 
marizes the harmonic oscillator quantum numbers for the pair 
coordinate. The oscillator parameters for the relative and the 
pair motion are given by 

^rel = 3/2 b, f>pair = ^ 2(A -To) ^ 


respectively. 

A compact notation can be obtained by defining the two- 
body density operator, which acts in two-body space, as 

R«= 2 \qQ)p a qQ , q , & (q'Q'\. (7) 

qQ,q’Q' 

The expectation value of any ‘bare’ two-body operator B 
can then simply evaluated as 

B a := = Tr 2 (r o B) , 


( 8 ) 
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where Tr 2 denotes the trace in two-body space. If one evolves 
the observable B in the same way as the Hamiltonian, namely 
B a = UlBU a , the expectation value 

0F a |B ff |'F o > = 01 • a \uiBU a \'¥ a ) 

= <'F a=0 |S|Y a=0 ) 

does not depend on a because = |'P a . = o). However, as 

for the Hamiltonian, the evolved observable B a is in general 
no longer a two-body operator and contains induced higher- 
body operators. 

The two-body approximation consists in calculating 

B a := Tr 2 (r o B a ) = Tr 2 (uJjJj r B) , (10) 

where the two-body operator B a is SRG transformed in two- 
body space. It should be noted that U a acts only on the rela¬ 
tive coordinate part of R a or B and leaves the center of mass 
motion of the pair unchanged. 

If the two-body approximation was exact (for both the 
Hamiltonian H and the observables B ) B a would be a- 
independent. The remaining a dependence therefore indicates 
the size of the neglected contributions from the induced three- 
and higher-body terms. 

With the many-body states V F„) obtained in the NCSM var¬ 
ious bare and SRG transformed two-body quantities will be 
studied in the following. With bare operators we can investi¬ 
gate how the properties of the eigenstates |'F ( ,) change with 
increasing flow parameter. This will reflect the increasing 
‘softness’ of the transformed Hamiltonian II On the other 
hand properties calculated with transformed operators should 
be independent from the flow parameter if the two-body ap¬ 
proximation is justified. However one has to be careful when 
drawing conclusions about omitted higher-order terms in the 
transformed operators. In this paper the two-body approxima¬ 
tion is employed both for the transformed Hamiltonian and for 
the transformed operators. Even if the two-body approxima¬ 
tion is perfect for the transformed operators we would still ex¬ 
pect a dependence on the flow parameter due to the two-body 
approximation for the transformation of the Hamiltonian. 

It might be possible however that particular components 
of the wave function are less sensitive to higher-order terms 
in the transformed Hamiltonian and can be obtained reliably 
within the two-body approximation. This appears to be the 
case for the pairs with small pair momentum as will be dis¬ 
cussed in Sec. HIE. 



a [fm 4 ] 

FIG. 1. (Color online) Contributions from the different S, T channels 
to the total energy of 4 He as a function of the flow parameter a for 
the AV8’ interaction (solid) and the N3LO interaction (dashed). The 
total energy E tot (black lines) includes the Coulomb energy. 


three- and four-body interactions (see Eq. (3)) that are ne¬ 
glected when calculating the 4 He ground state energy in two- 
body approximation. As seen in Fig. 1 the total energy varies 
by about 10% indicating that the contribution of three- and 
four-body terms are small. In Fig. 1 we also show the indi¬ 
vidual contributions from the four S,T channels to the total 
ground state energy of 4 He. These can be calculated by using 
Eq. (10) 

E a , ST = Tr 2 (/G[(f int + yf ] )n sr ®7]) , (11) 

where the operator n. 57 - projects on the spin-isospin channel 
S, T. The total energies for the AV8’ and N3LO interactions 
are very similar (although the individual kinetic and poten¬ 
tial contributions are quite different) and in both cases the 
dominant contribution to the binding energy is coming from 
S,T = 1,0 pairs with a large contribution from the tensor 
force. The S, T = 0,1 channel gives a much smaller attrac¬ 
tive contribution whereas the S,T = 1,1 channel provides 
repulsion and the contribution from the S,T = 0,0 channel is 
negligible. With increasing flow parameter we observe a re¬ 
duction of the attractive contribution of the S, T = 1,0 chan¬ 
nel and the S,T = 1,1 contribution becomes less repulsive. 
These changes are related to changes in the occupation of the 
different S, T channels that will be discussed in Sec. HID. 


B. Relative density distribution in coordinate space 


III. RESULTS 
A. Energies 

In two-body approximation U a is exactly unitary in two- 
body space, but in many-body space unitarity is only approx¬ 
imate. Therefore the energy eigenvalue E a depends on a and 
this dependence can be taken as a measure for the induced 


The repulsive nature of the nuclear interaction is easily seen 
in the the relative density distribution p\f(r) that gives the 
probability to find a pair of nucleons at a distance r. It is cal¬ 
culated by replacing B in Eq. (8) with the two-body operator 


V \rlm,SM s ,TM T )(rIm,SM s ,TM T \ ® 1 , (12) 

ImS Ms T Mj 
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FIG. 2. AV8’ (left) and N3LO (right) two-body densities in coordi¬ 
nate space calculated with bare (top) and SRG transformed density 
operators (bottom), a in units of fm 4 . See also Fig. 1 in Ref. [43] for 
the AV8’ two-body densities in the different S, T channels. 


where we sum over all quantum numbers except the radial 
distance r. 

The short-range correlations induced by the nuclear inter¬ 
action are reflected in the relative density distributions p]f(r) 
shown in the upper part of Fig. 2. They are calculated with 
the eigenstates of the bare and SRG evolved Hamiltonians H a 
according to Eq. (8). For the bare interactions the relative den¬ 
sity distributions p ru ' (r) show the typical correlation hole at 
short distances. Due to the repulsive core of the interaction 
the probability to find a pair of nucleons at a distance r be¬ 
tween their centers is depleted at distances below r as 1 fm. 
The chiral N3LO interaction is obviously not as repulsive as 
the AV8’ interaction at short distances, as the correlation hole 
is less pronounced for the N3LO interaction. 

One should keep in mind that the term ‘correlation hole’ 
is misleading as for r ~ 0.5 fm the nucleon densities are al¬ 
ready strongly overlapping so that there is no hole in the bary- 
onic matter density. On the contrary at short distances one 
encounters locally large baryonic matter densities and expects 
strongly polarized nucleons. 

With increasing flow parameter a the correlation hole dis¬ 
appears more and more. Furthermore the density distributions 
calculated for the evolved AV8’ and N3LO interactions be¬ 
come increasingly similar, cf. Fig. 2 (a), (b). For the largest 
flow parameter the relative density distribution is essentially 
of Gaussian shape — as would be expected for an uncorre¬ 
lated mean-field wave function. Thus, the SRG transforma¬ 
tion brings us from a highly correlated system to a simple shell 
model or mean-field like situation. 

On the other hand the density distributions p^'fr) shown 


in the lower part of Fig. 2, which are obtained with the 
SRG transformed density operators according to Eq. (10), are 
all very similar to those obtained for the bare Hamiltonians 
H = H a=0 . Without employing the two-body approximation 
U a would be unitary in four-body space and p a would not de¬ 
pend on a. The remaining a dependence indicates that the in¬ 
duced three- and four-body terms are small but not completely 
negligible. 

The probability densities obtained with both bare and trans¬ 
formed density operators are normalized to the number of 
pairs: 

^ drr 1 pf(r) = £ dr r 2 p?(r) = A(A ~ ' j (13) 

The spatial distributions shown in Fig. 2 may help our in¬ 
tuition but can not easily be related to experiment. Therefore 
momentum distributions and momentum correlations will be 
addressed in the following. 

C. Relative momentum distributions 

The relative momentum distribution, i.e., the probability 
n^i ST (k) to find a nucleon pair with relative momentum k, rel¬ 
ative orbital angular momentum /, spin S, and isospin T is 
obtained by replacing B in Eq. (8) with the two-body operator 

^r(k) = 

Y \klm,SMs, TM T )(klm,SMs, TM t \ ®T, (14) 

mMs Mj 

where \klm,SMs, TMj) denotes the spherical momentum 
space representation of the relative motion. 

The relative momentum distributions for pairs with relative 
momentum k, irrespective of their orbital angular momentum 
l, spin S and isospin T, n r c f{k.) = YjIst n a\sr(k), and its SRG 
transformed partner h l ^{k) are shown in Fig. 3 for the two 
interactions and different flow parameters. The short-range 
repulsive correlations, which manifest themselves in coordi¬ 
nate space as the correlation holes, show up as tails in the 
momentum distributions n^ l (k) that reach out to large relative 
momenta k. Note however that the momentum distribution is 
not the Fourier transform of the diagonal two-body density in 
coordinate space. 

Whereas the tail of the momentum distribution shows an 
exponential behavior at large relative momenta for the AV8’ 
interaction, the N3LO relative momentum distribution reflects 
the momentum space regulator that cuts off high momenta be¬ 
yond about 3.5 for 1 . For both interactions short-range ten¬ 
sor correlations play an important role as will be discussed in 
more detail later. With increasing flow parameter a the high- 
momentum components are more and more reduced until the 
probability distribution of the relative momentum assumes a 
Gaussian shape for both evolved interactions, cf. Fig. 3 (a) 
and (b), corresponding to an uncorrelated wave function. 

Fig. 3 (c) and (d) show that the density distributions h^ x (k) 
obtained with the SRG transformed density operators are 
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FIG. 3. AV8’ (left) and N3LO (right) two-body densities in 
momentum-space calculated with bare (top) and SRG transformed 
density operators (bottom), a in units of fm 4 . See also Fig. 2 in 
Ref. [43] for the AV8’ two-body densities in the different S, T chan¬ 
nels. 


again all very similar to those obtained for the bare Hamil¬ 
tonian, indicating that the induced three- and four-body terms 
are also small in momentum space. The most visible depen¬ 
dence on a occurs for the AV8’ interaction around k ~ 2 fm 1 . 
This dependence is related to three-body correlations induced 
by the two-body tensor force to be discussed in the following 
subsection. 


D. S,T relative momentum distributions 

In order to get a deeper understanding of the nature of the 
correlations we separate the momentum distribution accord¬ 
ing to Eq. (14) into its parts h^ ST (k) = Z/ com i n g 

from the different S,T channels and display the results in 
Fig. 4. The first observation is that the S,T = 1,0 channel, 
in which the tensor force is strongest, shows only a weak a 
dependence. In contrast to that there is strong a dependence 
in the S , T = 0, 1 channel for momenta around 2 fm 1 and 
also in the S, T = 1,1 channel. With increasing a strength 
moves from the odd S, T = 1,1 to the even S , T - 0, 1 chan¬ 
nel. This effect has been discussed in detail in Ref. [37]. In 
a simplified picture one may consider the situation of a local¬ 
ized S, T = 0, 1 pair with a third nucleon of different flavor 
not too far away. Additional binding from the strong tensor 
force acting between the third nucleon and a nucleon of the 
pair can be obtained by flipping the spin of one nucleon in the 
S, T = 0,1 pair and thus converting with some probability the 
S, T = 0,1 pair into an S , T = 1,1 pair. It is energetically 
more favorable to break some S, T = 0,1 pairs and loose their 


TABLE I. Number of pairs in the different (5, T) channels for the 
bare and SRG evolved AV8’ and N3LO interactions. 


pair ~pair 

n u$T ~ n a,ST 

(0,0) 

(0,1) 

(1,0) 

(U) 

AV8’, bare 

0.008 

2.572 

2.992 

0.428 

AV8’, a = 0.01 fm 4 

0.008 

2.708 

2.992 

0.292 

AV8’, a = 0.04 fm 4 

0.007 

2.821 

2.993 

0.179 

AV8’, a = 0.20 fm 4 

0.005 

2.925 

2.995 

0.075 

N3LO, bare 

0.009 

2.710 

2.991 

0.290 

N3LO, a = 0.01 fm 4 

0.007 

2.745 

2.992 

0.255 

N3LO, a = 0.04 fm 4 

0.006 

2.817 

2.994 

0.183 

N3LO, a = 0.20 fm 4 

0.004 

2.921 

2.995 

0.079 


binding if more binding from the tensor interaction with the 
third particle can be gained. These genuine three-body corre¬ 
lations are lost in two-body approximation and introduce an 
a-dependence. With increasing a the tensor part of the SRG 
transformed Hamiltonian H a is weakened while the central 
part is strengthened. With the weakening of the tensor force 
the effect of the three-body correlations will be reduced. This 
explains the reduction of the number of S , T = 1,1 pairs with 
increasing a and the a-dependence of the S,T = 0,1 pairs at 
momenta around 2 fm 1 where the relative importance of the 
tensor force is most pronounced. 

This transfer of probability between the even S, T = 0,1 
and the odd S,T = 1,1 channel can be seen in Table I where 
the number of pairs, n^ T given by the integrals over the mo¬ 
mentum distributions, are listed as a function of a. With in¬ 
creasing flow parameter the occupation numbers approach the 
limit of the mean-field or independent particle model with 3 
pairs in both even channels and 0 pairs in the odd channels. 
The number of pairs n^ T and n^ T obtained with bare and 

SRG transformed density operators are identical as U a in two- 
body approximation does not connect different S, T channels. 

The reshuffling of probability between the different spin- 
isospin channels also tells us that the omitted many-body 
terms in U a and the transformed Hamiltonian //„ will have 
a non-trivial spin- and isospin-dependence. 


E. Relative momentum distributions for K = 0 pairs 

Up to now we have investigated relative momentum distri¬ 
butions for all pairs, indiscriminate of the pair momentum K. 
It has been found that the relative momentum distributions de¬ 
pend quite significantly on the pair momentum [23, 44]. In the 
context of this paper it is interesting to see how this is related 
to many-body correlations. 

We might expect that back-to-back pairs with K « 0 are 
less affected by many-body correlations than pairs with a large 
pair momentum K. In a K = 0 pair with large relative momen¬ 
tum k both nucleons have large individual momenta. For pairs 
with large pair momentum K however there is a high proba¬ 
bility that one of the nucleons has a momentum less than or 
close to Fermi momentum. We would therefore expect that 
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FIG. 4. (Color online) AV8’ (top) and N3LO (bottom) relative momentum distributions in the different spin-isospin channels obtained with 
SRG transformed density operators, a in units of fm 4 . 


these nucleons are interacting more strongly with other nucle¬ 
ons and therefore are susceptible to many-body correlations. 

In order to study this we investigate the more exclusive joint 
probability h a ,sT(k, K) to find a nucleon pair with spin S and 
isospin T at relative momentum k and total pair momentum 
K. It is calculated with the two-body operator 

nisr(k, K ) = 

^ \klm, S M s , TM r ){klm, S M s , TM r \ 

mMs Mj 

®Y} KLM ^{j^2j V2 ( klm \’ (! 5 ) 

LM 

where \KLM) denotes the spherical momentum representa¬ 
tion of the relative momentum of the pair with respect to the 

/ , \3/2 

A - 2 remaining nucleons. The factor I J originates from 
the transformation from Jacobi coordinates to the coordinates 
k and K, see Sec. D. In this paper we consider only pair 
momentum K = 0 for which the sum over L, M reduces to 
L = 0,M = 0. 

The momentum distributions n a £T(k,K = 0), obtained by 
summing h a jsT(k,K = 0) over /, are displayed in Fig. 5. One 
sees that in all S, T channels the results are essentially inde¬ 
pendent on a. We can also observe that with increasing flow 
parameter a the momentum distributions ii a ,ST(k), as shown 
in Fig. 4, become more and more similar to the momentum 
distributions n a sT(k,K = 0) . This consistently tells us that 
for pairs with total momentum K — 0 many-body correlations 
are not very important, as anticipated in the discussion above. 
The K = 0 pairs are therefore the best candidates for experi¬ 
mental studies of short-range two-body correlations. Similar 
considerations can be found in Refs. [12, 23, 44]. 


We can also notice significant differences between the two 
even channels. In the S, T - 0, 1 channel the momentum dis¬ 
tribution has a node at relative momenta of about 1.8 far 1 . 
This is very different in the S, T = 1,0 channel. Here the mo¬ 
mentum distribution does not show a minimum and the num¬ 
ber of pairs for relative momenta above about 1.5 far 1 is sig¬ 
nificantly larger. This difference is due to short-range tensor 
correlations that only contribute in the S,T = 1,0 channel. 
To illustrate this we show in Fig. 6 the momentum distribu¬ 
tions in the S , T — 1,0 channel decomposed into their contri¬ 
butions from different / for the AV8’ and N3LO interactions. 
For the K — 0 distributions we do not have the results for the 
bare AV8’ interaction and used the SRG transformed densi¬ 
ties n r ^ T (k, K = 0) for a = 0.01 fm 1 instead. Because of 
the very weak a dependence this should be equivalent to the 
exact result for the bare AV8’ interaction. The distributions 
n aJST (k,K = 0) in this channel are of particular interest, be¬ 
cause they show the dominant contribution from pairs with 
relative angular momentum l = 2 in the momentum region 
from 1.5 to almost 4 fm for the AV8’ and from 1.5 to about 
2.5 far 1 for the N3LO interaction. These D -wave pairs are di¬ 
rectly reflecting the correlations induced by the tensor force. 
In the AV8’ case the tensor correlations are present at all rel¬ 
ative momenta and dominate over the / = 0 contribution for 
momenta above k ~ 1.5 fm 1 . For N3LO the regularization 
cuts them off at higher momenta so that they dominate only 
between^ ~ 1.5 far 1 and A: ~ 2.5 far 1 . It is important to note 
that the tensor correlations contribute also at low momenta. 
There is no natural scale in the tensor correlations that would 
allow to separate low and high momentum regions. 

It is also interesting to compare the 4 He momentum distri¬ 
butions in the S,T - 1,0 channel with those of the deuteron 
























































7 


n 


10 3 

AV8' s-o,t-o 

10 2 

- - (*=0.01 . 


- (*=0.04 

10 1 

-— ar=0.20 ' 

10° 


1CT 1 



f \ 

1CT 2 

1 \ 

1CT 3 

(a) \ 


\ 


\ 


\ 


S=0,T=1 

- - or=0.01 

- (*=0.04 

■-(*= 0.20 


(f» 


\c 



S=1 ,T=0 i 

/ 

0 

- - (*=o.oi ; 

- (*=0.04 

- (*=0.20 1 

.(c) 




S=1 ,T=1 
- - (*=0.01 . 

- (*=0.04 

-— (*=0.20 ' 

r\ 

\\ 

( d) w 




FIG. 5. (Color online) AV8’ (top) and N3LO (bottom) relative momentum distributions for pairs with vanishing pair momentum K = 0 in the 
different spin-isospin channels obtained with SRG transformed density operators. 


shown in Fig. 7. Whereas the momentum distributions in the 
low momentum region up to about 1.5 fm 1 are noticeably 
different, reflecting the differences in the long-range parts of 
the wave functions for the loosely bound deuteron and the 
strongly bound 4 He, the momentum distributions are almost 
indistinguishable for momenta above 1.5 far 1 . This is consis¬ 
tent with the observed universality of short-range correlations 
discussed in Ref. [37]. 

If we compare the momentum distributions for the S, T = 
1,0 pairs integrated over all pair momenta K as shown in 
Fig. 6 (a) and (b) with those for the corresponding K — 0 
pairs we notice that the node for the S -wave pairs has van¬ 
ished and that we find additional contributions from higher 
relative orbital angular momenta (G-wave pairs). This is con¬ 
sistent with our expectation that many-body correlations play 
a greater role for pairs with large pair momentum K. This is 
further confirmed by a strong a dependence for these momen¬ 
tum distributions (not shown here). Alvioli et al. [44] found 
that the AV8’ momentum distributions factorize in a relative 
and a pair distribution for pair momenta K up to about 1 far 1 , 
whereas for higher pair momenta the momentum distributions 
nsr( k, K) no longer factorize and depend on the relative ori¬ 
entation of k and K. 


F. Relative probabilities for S, T pairs 

The strength of two- and many-body correlations strongly 
depends on the relative momenta of the pairs. To highlight the 
relative importance of these correlations it is advantageous to 


look at the relative probabilities 

n’sj-(k) n ST (k,K - 0) 

lsmf T (k) ’ 1st n ST (k, K = 0) 

to find pairs in a given S, T channel as a function of relative 
momentum k. In Fig. 8 we show these relative probabilities 
for the bare AV8’ and N3LO interactions. If we look at the rel¬ 
ative probabilities for all pairs with relative momentum k it is 
obvious that the S, T = 1,1 pairs contribute significantly, es¬ 
pecially around k ~ 1.5 far 1 . The total number of S,T = 1,1 
pairs is much smaller than the number of pairs in the even 
channels as shown in Tab. I, but in this mid-momentum re¬ 
gion the number of pairs in the S, T = 1,1 channel is compa¬ 
rable to those in the even channels. This is the region where 
many-body correlations have the largest effect. For small rel¬ 
ative momenta the relative probabilities are dominated by the 
mean-field, for relative momenta above about 3 fur 1 in case 
of the AV8’ and about 2.5 fm 1 in case of the N3LO inter¬ 
action the relative probabilities are dominated by short-range 
correlations. This influence of many-body correlations is also 
related to a strong dependence on a. With increasing flow pa¬ 
rameter a the relative probabilities for all pairs become more 
and more similar to the relative probabilities of the K — 0 
pairs shown in Fig. 8 (c) and (d). As there is no significant a 
dependence for the K = 0 momentum distributions, the rela¬ 
tive probabilities for the K — 0 pairs are independent from a 
as well and are therefore not sensitive to many-body correla¬ 
tions. 

It is interesting to note that the relative probabilities for 
K - 0 the pairs are quite similar for the AV8’ and N3LO in¬ 
teractions, even for very large relative momenta, whereas the 
absolute values of the momentum distributions are very dif- 
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k [fm -1 ] 


k [fm -1 ] 


k [fm -1 ] 


k [fm -1 ] 


FIG. 6. (Color online) Momentum distributions for the bare AV8' 
and N3LO interactions in the S,T = 1,0 channel decomposed into 
contributions from pairs with relative orbital angular momentum 
1 = 0,2,4. For all pairs (top) and for pairs with vanishing pair mo¬ 
mentum (bottom). 



FIG. 7. (Color online) Total, / = 0 and 1 = 2 momentum distributions 
in the deuteron (S, T = 1,0) for the AV8’ and the N3LO interactions. 


ferent. Differences in the relative probabilities between AV8’ 
and N3LO reflect differences in the relative importance of ten¬ 
sor correlations for the two interactions due to differences in 
the regularization of the tensor force. 


G. Relative probabilities for pn and pp pairs 


In an experiment one measures protons and neutrons and 
not S , T pairs. Therefore we define the operator in two-body 
space that measures the probability to find a pair of two pro¬ 
tons (pp: Mj = 1) or a proton-neutron pair (pn + np: Mj = 0) 


FIG. 8. (Color online) Relative probability to find pairs in the differ¬ 
ent spin-isospin channels as a function of relative momentum. For 
all pairs (top) or only for pairs with pair momentum K = 0 (bottom). 


at relative momentum k and pair momentum K as 


n Mr (k, K) = 

| klm,SM s , TM T )(klm, S M s , TM t \ 


lm,SM s ,T 


/ A \ 3/2 
Y)KLM ){—) (KLM |. (17) 


In the case of 4 He one T = 0 pair corresponds to one pn- 
pair, and one T — 1 pair to | of a pp, | of a tin and | of a pn 
pair. 

The corresponding relative probabilities for K = 0 pairs 
n Mr (k, K — 0) 

-( 18 ) 

YjM t nM T (k,K = 0 ) 

are shown in Fig. 9. The first observation is that the relative 
probabilities are rather similar for the AV8’ and N3LO inter¬ 
actions. At low momenta both show a ratio close to | to find 
pp versus pn pairs. This is to be expected because an uncor¬ 
related system of two protons and two neutrons can form one 
pp- pair, one nn-pair but four pn-pairs. 

Around k ® 1.8 fm -1 the minimum in the S , T = 0,1 chan¬ 
nel (see Fig. 5) together with the / = 2 contribution from 
the tensor interaction in the S,T = 1,0 channel (see Fig. 6) 
enhances the relative probability to find a pn pair to almost 
100%. This dominance of pn pairs has been observed in ex¬ 
clusive two-nucleon knockout experiments [13, 15]. Recently 
the pn to pp ratio has been measured for relative momenta k 
from 2.5 far 1 up to almost 4 far 1 [17] showing an increase 
in the pp/pn ratio. Within the experimental uncertainties the 
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FIG. 9. (Color online) Relative probability to find pn or pp pairs 
with pair momentum K = 0 as a function of relative momentum. 


data agree with our results for both AV8’ and N3LO interac¬ 
tions. 


IV. SUMMARY AND CONCLUSIONS 

In this paper we applied the SRG formalism for the calcu¬ 
lation of relative density and momentum distributions of 4 He. 
The 4 He ground state wave functions are calculated in the 
NCSM with the SRG evolved AV8’ and N3LO interactions 
in two-body approximation. Two-body densities in coordi¬ 
nate and momentum space calculated with the unevolved den¬ 
sity operators illustrate how short-range correlations are elimi¬ 
nated by the SRG evolution. With increasing flow parameter a 
the interaction gets ‘softer’ and the wave functions become es¬ 
sentially uncorrelated mean-field wave functions without cor¬ 
relation holes and high-momentum components. The short- 
range or high-momentum information can be recovered by 
calculating two-body densities with the SRG evolved density 
operators, again in two-body approximation. Using these ef¬ 
fective density operators we see a dependence of the two-body 
densities on the flow parameter a. This a dependence is due 
to missing contributions from three- and four-body terms in 
the effective operators, that are omitted in the two-body ap¬ 
proximation. We find that not all components of the two-body 
density are equally affected by many-body correlations. The 
momentum distributions for pairs with pair momentum K = 0 
show only a very weak a dependence and therefore provide 
direct access to two-body short-range correlations. The mo¬ 
mentum distributions integrated over all pair momenta on the 
other hand have a sizeable a dependence. This a dependence 
is particularly strong in the S,T — 0,1 and S,T = 1,1 chan¬ 
nels for momenta around 1.5 fin -1 . We identified three-body 
correlations induced by the strong tensor force as the main 
contributor to these many-body correlations. 

Our results for the AV8’ momentum distributions agree 
with previous results [22-24, 44] using variational Monte- 
Carlo and few-body approaches. We confirm the important 
role of tensor correlations that explain the experimentally 
observed dominance of pn over pp pairs in exclusive two- 
nucleon knock-out with large momentum transfer. We also 
show that the chiral N3LO interaction provides very similar 


results for the momentum distributions at least up to relative 
momenta of 2.5 - 3.0 fm -1 . This includes the role of tensor 
correlations that are very similar for both interactions. At 
larger relative momenta the N3LO momentum distributions 
fall off much faster than the AV8’ momentum distributions as 
would be expected from the relatively soft cutoff employed 
in the regularization of the N3LO interaction. These differ¬ 
ences are however mostly hidden in the ratios of pn versus 
pp pairs as a function of relative momenta as investigated 
in two-nucleon knockout experiments. Recently new chiral 
interactions with different regularization schemes have been 
proposed [45, 46]. It will be interesting to see how these af¬ 
fect the short-range behavior and momentum distributions. 

In this paper we investigated the 4 He nucleus as it allowed 
us to compare SRG evolved with bare interactions. SRG 
evolved soft Hamiltonians however will allow us to study 
heavier nuclei in the p-shell using the NCSM. One interesting 
question is the isospin dependence of the short-range correla¬ 
tions [47, 48] in asymmetric nuclei. Neutron halos or skins 
might be a useful laboratory for the study of neutron mat¬ 
ter. The present study also did not include three-body forces. 
Wiringa et al. [24] find only small differences in the proton 
momentum distributions obtained with AVI8 alone and AVI8 
together with an UX three-body interaction. We would expect 
that three-body interactions would not significantly change the 
short-range two-body correlations in the K = 0 pairs, they 
may however have a significant effect on the many-body cor¬ 
relations and therefore change the results for the momentum 
distributions integrated over all pair momenta. 


Appendix A: Talmi-Moshinsky transformation 


For an orthogonal transformation of the coordinates (xi, X 2 ) 
into (X, x) with the mass ratio d = A 1 /A 2 


X' 

.x. 




(Al) 


the product of harmonic oscillator wave functions in coordi¬ 
nates Xi and X 2 can be expressed with the Talmi-Moshinsky 
brackets (NL,nl: A| «i/i, 712/2 ■ A) d in the new coordinates 
as: 


z 


h h 

m\ m 2 


^ < NL,nl : A| nili,n 2 h : A ) d 

z 


NL,nl 


Mm 


L l 
M m 


^lm(X)^1(x) (A2) 


All harmonic oscillator wave functions have the same os¬ 
cillator parameter b. See Ref. [49] for further properties. 
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Appendix B: Jacobi coordinates 


and the conjugate coordinates in momentum space as 


The wave function and two-body densities in the intrinsic 
system will be expressed in Jacobi coordinates. We follow 
[50] and use mass-scaled Jacobi coordinates that fulfill the or¬ 
thogonality condition (Al): 


£o - 
€1 = 
£2 = 


- [Xj + . . . + X A ] 

A 




■X 2 ] 


^(Xi +x 2 )-x 3 


(Bl) 

(B2) 

(B3) 


k =^<Pi-P2)> K = p,+p 2 . (C4) 

These coordinates differ from the ones in (Al) by factors of 
V 2 and 1 / V 2 so that the oscillator parameters for the relative 
and the pair motion have to changed to b K \ = V 2 b and b pa ; r = 
-q= b respectively. Using (A2) and coupling orbital angular 
momenta, spins and isospins the two-body density (C2) can 
be transformed into P^q ^q,- 
Using this new basis 

\qQ) = l^rei; nlmSM s TM t ) ® \b piur -NLM) (C5) 


we can define the density operator in two-body space as 


£ 4-1 


lA- 1 

1 , 

/ A 

A j(Xl + ... +X A _!) X A 


(B4) 


2 \qQ)p l X- & WQ'\- (C6) 

qQ-q'Q' 


The transformation from the coordinates (xi^.-.x^) to 
(g 0 ,... ,£a-i) ' s orthogonal. Translational invariant one-body 
densities can be expressed in these Jacobi coordinates. For the 
two-body densities we prefer however a different set of Jacobi 
coordinates where we have the distance between the nucleons 
and the distance of the center of mass of the pair with respect 
to the rest of the nucleus as coordinates. This can be achieved 
by an additional orthogonal transformation from the coordi¬ 
nates (£ 4-21 £ 4 - 1 ) > nt0 the coordinates (//, it) 


/? ,lh allows to easily express the two-body density in coordi¬ 
nate space as a function of relative distance r and pair position 

R 

p lab (r, R) = Tr 2 (fl* ab |r; R)<r; R|) , (C7) 

or as function of relative momentum k and pair momentum K 
7 ? lab (k, K) = Tr 2 (^ ab |k; K)<k; K|) . (C 8 ) 


2 (A - 2) 


[^-^(xi + ... + x A _ 2 ) - ^(xa-i + x A )], (B5) 





(B 6 ) 


Appendix C: Two-body densities in the laboratory system 


This can be trivially extended to calculate off-diagonal densi¬ 
ties or densities for pairs of a given spin and isospin. 

These densities can also be expressed as 

p lab (r, R) = 

-xj - r)tf 3 (^& + xj) - R)|'T) (C9) 

i<j 

and can be obtained from the wave functions by integrating 
out A - 2 coordinates (we have omitted spin- and isospin in¬ 
dices for brevity) 


In the NCSM the wave function is expanded in a harmonic 
oscillator single-particle basis with oscillator parameter b 

| q) = \b',nlm,m s ,m t ) . (Cl) 

The NCSM uses second quantization techniques and the 
two-body density in the harmonic oscillator basis can be ex¬ 
pressed as 

= WKKwl'*) ■ (C2) 

For the discussion of two-body correlations it is natural to 
transform from single-particle to pair coordinates. In the lab¬ 
oratory system we define relative and pair coordinates as 

r = X!-x 2 , R=^(xi+x 2 ), (C3) 


p lab (r, R) = A( " A ^ — J" d 3 xi ■ ■ ■ (Pxa -2 

'J , (x 1 ,...,x A _ 2 ,R+ir,R-ir)‘ 

x¥(xi. x a . 2 ,R+ 3 UR- ±r), (CIO) 

where we have used the antisymmetry of the wave function. 

Appendix D: Two-body densities in the intrinsic system 

The two-body density p' a g ,g, is calculated in the harmonic 
oscillator basis of the NCSM localized at the origin of the co¬ 
ordinate system and is therefore not translationally invariant. 
However the wave function in the NCSM factorizes into an in¬ 
trinsic wave function only depending on relative coordinates 
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and the total center-of-mass wave function in the ground state 
of the harmonic oscillator [5 1 ] 

nxu ... ,XA) = .. ,£ a _ 3 , j/,#)^oo(^o) • (Dl) 

This allows us to relate the two-body density in the labo¬ 


ratory system with the two-body density in the intrin¬ 

sic system. The derivation is lengthy but straightforward and 
follows the derivation of the translational invariant one-body 
density in Ref. [50]. 

One starts with the two-body density in the laboratory sys¬ 
tem as given from the wave functions written in the coordi¬ 
nates Xi,. .. ,x A _ 2 .X 2 = ^(x A _!+x A ) and # = -^(x A _!-x A ): 


lah 

Pq' ,N'L'M'-qN 2 L 2 M 2 


A(A - 1) 


f 


d 3 x\ 


■ d 3 x A -2 d 3 X 2 d 3 d 


f d>x 'i 


d 3 x A ~ 2 d 3 X' 0 d 3 # 


x 'T(xi,... ,x a _ 2 ,X 2 ,i?)*(5 3 (xi - Xj) • • • c5 3 (x a ^> - x' A _ 1 )'V{x\,... ,x' A _ 2 ,X 2 ,#) 


x ^V'Z.'M;(X 2 )* <PN 2 L 2 M 2 (X 2 ) <fq’(#) ‘Pqift) (D2) 


We now perform an orthogonal coordinate transformation from (xi,..., x A _ 2 , X 2 ) to (X A _ 2 , ,..., £ 4 . 3 ) and correspondingly 

for the primed coordinates. Next one uses the properties of the delta function to obtain delta functions in the Jacobi coordinates 
ij k and in X A _ 2 = —== Ef ] 2 x,. The delta function in X A 2 is expanded in the harmonic oscillator basis 


£ 3 (X A _ 2 


X A - 


») - X ‘PNa-iLa- 2 M a _2 (X A - 2 ) <PNa- 2 .La- 2 .Ma -2 (X a _ 2 )* 


Na—2 La~2 Ma~~> 


(D3) 


In the next step a second orthogonal coordinate transformation from (X A _ 2 , X 2 ) to (£ 0 , 17 ) is performed employing the Talmi- 
Moshinsky transformation for the harmonic oscillator wave functions, rewriting products iPn a - 1 La- 2 M a ^('X-a- 2 )‘Pn 2 L 2 M 1 (X- 2 ) with 
linear combinations of products <PN cm L an M cm (^o)‘PNLM(ii)- 

If one uses (Dl) one can now integrate over £ 0 and €0 t0 express the density matrix in the laboratory system with the intrinsic 
wave functions so that the density matrix in the laboratory system can be related to the density matrix in the intrinsic system: 


Pq',N'L'M';qNLM ~ 


A(A - 1) 


x 'T inl (^|,..., f A _ 3 , 1 h &)*% .., ^f A _ 3 , if,#) <Pn.um.WT <PnlmW) <p q &) 0>4) 


In this paper we only discuss scalar densities in the pair coordinates (integrating over all pair momenta K or K = 0). Then 
L — L' and M — M' and we can use the completeness relations for the Clebsch-Gordan coefficients to obtain the relation between 
the density matrices in the laboratory (D2) and the intrinsic system (D4) 


V> b 

/ , Pq'N' 2 L 1 Mz,qN 7 .L 1 M 1 
M 2 


-E 

NN'L 


M 


( 0 ) 


Z< 

M 


NiLi.NiL->:N'L.NL 7 ,Pq'N'LM-,qNLM 


(D5) 


with the transformation matrix M (()) 

M nLn 2 l 2 -,n'l,nl = X < 00 ’ NL ■ L| Na-2L A -2,N 2 L 2 : L) (A _ 2)/2 <00, N'L : L\ N A - 2 L A - 2 , N' 2 L 2 : L) (A _ 2y2 (D 6 ) 

Na-iLa-2 

given by the Talmi-Moshinsky brackets with the mass ratio d = This reflects the distribution of the oscillator quanta among 
the nucleon pair and the A-2 remaining nucleons. The translationally invariant two-body density averaged over the orientations 
of the pair momenta Em Pq'N'LM qNLM is then obtained from the two-body density matrix in the laboratory system by inverting 
Eq. (D5). 

The two-body density in the Jacobi coordinates j] and 1 ? is given by 

//) = X Pq'N'L'M'-qNLM <PnlmW) Vn'L'M'W) <Pq($)* Vq'i#) (D7) 

qNLM.q'N'L'M' 


The mass-scaled Jacobi coordinates have technical advantages, however we prefer to express the two-body densities in the 
more intuitive coordinates 


r = x A _!-x A , R = -(x A _! +x A ) - j—- 


(Xi + . . . + X A _ 2 ) 


(D 8 ) 
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that are related to the Jacobi coordinates by 


r = V2tf, R = 


2 (A - 2) 




(D9) 


The conjugate variables to r and R are the relative momentum of the nucleons in the pair k and the relative momentum of the 
pair with respect to the rest of the nucleus K (in the center-of-mass system this is the same as the momentum of the pair). 


1 A-2 2 

k = 2 (Pa-i - Pa), K = —(p A _r + p A ) - ^(Pi + ■ ■ ■ + p A _ 2 ) 


(DIO) 


The two-body density in these coordinates is then 




A 

A- 2 


3/2 


V2 

z 

qNLM,q'N'L'M' 


Pq'N'L?M'\ q NLM ^fc M ,(R) ^ l 2 3 4 5 6 7 8 9 10 (r) 


(Dll) 


Note that the oscillator parameter for the pair motion b is different than in the laboratory system, whereas the 

oscillator parameter b le \ = V2 b is the same. The factor ^ is the Jacobian of the nonorthogonal transformation from the 
coordinates j'J and // to r and R in (D9). 


As in the laboratory system a convenient notation is ob¬ 
tained by defining the density operator in two-body space 

R = J] \<lQ)PqQ-,q'Q'WQ'\ ( D12 ) 

qQ.q'Q' 


and in momentum space as 


n{ k,K) = Tr 2 



with 


| qQ) = \b lel -,nlmSM s TM T ) ® |fc pair ; NLM) . (D13) 
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